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1. INTRODUCTION 

Heavy quark effective theory is routinely used in phenomenology. In these applications, 
the matching to QCD is achieved perturbatively and matrix elements of the operators in 
the effective theory are determined from experiment and models. However, HQET took 
its origin as an effective theory in the lattice regularization, where it was designed as a so- 
lution to the problem of treating quarks which are heavy compared to the inverse lattice 
spacing and thus do not propagate properly in the standard relativistic framework [H]]. 

Unfortunately, after considerable initial activity (see e.g. [0, S 0] and references 
therein) the non-perturbative treatment of the effective theory on the lattice had been 
somewhat dormant for a while. The reason is that it was realized [5] that a non- 
perturbative matching to QCD is needed; otherwise the continuum limit does not exist. 
A practicable solution of this problem was only found recently [6, 7, 8J. 

Here we point out that a non-perturbative matching is necessary on and ojfihe, lattice, 
the problem being most severe on the lattice and we review a non-perturbative test of 
HQET. We then explain a recent strategy to perform fully non-perturbative computations 
in HQET and discuss the status and perspectives of this approach. 
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2. HQET AS AN ASYMPTOTIC EXPANSION OF QCD 



Consider QCD at energies low enough such that the top-quark may be neglected alto- 
gether. In the QCD Lagrangian 

^QCD = -:^tr{F^yF^y}+Y,Wf[D^ir^i+mf]\j/f (1) 

^Sq f 

the sum over flavors then extends over / = u,d, s,c,b. An effective theory, HQET is 
expected to provide the asymptotic expansion of a certain (large) set of energies (e.g. 
mass splittings) and matrix elements of QCD in terms of the inverse of the mass of the 
b-quark.^ We restrict our discussion to the energies and matrix elements of states which 
contain a single b-quark at rest and refer to reviews such as [9, 10] for the general case of 
finite velocity. The Lagrangian of the theory, which may be obtained by a formal 1 /mb 
expansion (see e.g. 1 1 1]), is then given by the replacement 

Wh[D tirfi+mb]n ^stat + ^^^^ + --- , ^st^i = Wh[Do + 5m\\j/h. (2) 

Here ^^^^ is of order 1/mb and the mass term of the b-quark has been removed from 
the Lagrangian such that observable quantities in the b-sector have a finite limit as 
m\,^ oo (with a suitable counter term 5m). The effective heavy quark field i//h has only 
two degrees of freedom as appropriate for a non-relativistic spin 1/2 particle. Still it is 
notationally convenient to keep yji as a 4-component spinor but impose the constraint 

P+V/h = V^h, /'+ = (! + 7o)/2; (3) 

i.e. the lower components vanish in the Dirac representation. In order to discuss matrix 
elements, such as the B-meson decay constant, also the composite fields involving b- 
quarks are translated to the effective theory, for example: 

Ao = ZAV7iror5V^b ^ Ar = Zf^v7i7o75V/h. (4) 

Here Za,Z^^^ are the renormalization constants of the composite fields. 

The effective theory is valid for the low-lying energy levels as well as their matrix 
elements, the simplest one being 

4>QCD = ^ Za(0|Ao|5) . (5) 



It is scale independent, due to the chiral symmetry of QCD in the massless limit 
(including mt, = 0). In the effective theory this symmetry is absent and Zf^^ depends 
on the energy scale, /i, used in the renormalization condition which defines the finite 
current. Instead of 4>^'^'(jU) = Z^^'(/i) (0|Ao|5)stat it is therefore better to consider the 
renormalization group invariant matrix element 

4>RGi = lim [2b,g\^) ] 4>^'^V) ■ (6) 



^ Powers of 1 / mb are understood to be accompanied by slowly (logarithmically) varying functions of mb. 



It is both /i and renormalization scheme independent, as is easily seen using 
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by 



^'^u (u) 
scheme 



^s[,henie'(M)(l + In eq. ©, the coefficients bo.Yo defined 



P(s)=^^g = -bof + Oif) , rig) = 



II d 



^stat 



enter. We can now write down the HQET-expansion of the QCD matrix element 
$QCD = Cps(Mb/AMs)x4>RGi + 0(l/Mb), 
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Let us dicuss the somewhat unfamiliar form of eq. (HJ and the conversion function 
Cps(Mb/Apg). In a more conventional form we have 



4>QC° = CmatchK/M) x*Ms(M) + 0(l/mb) 



(11) 



with a matrix element 4>^(/i) renormalized in the effective theory in the MS scheme 
and the matching coefficient Cniatch('^b/M) depending on the b-quark mass mb in the 
MS scheme at scale mb, i.e. m^(mb) = mb. The factor Cmatch is determined (usually in 
perturbation theory) such that eq. (fTTT) holds for some particular matrix element of the 
current and will then be valid /or all matrix elements. Contact to eq. ^ is easily made 
by using 



^RGI 

%s(m) 
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MS 
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setting the arbitrary renormalization point /i to mb and identifying 



Mb 

^MS 



Cmatch ( 1 ; 
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^RGI 



(12) 



(13) 



= [^bQg (mb)] " "exp 



/match (g) To 



where g is taken in the MS scheme. The last equation may be taken as a definition of the 
anomalous dimension /match in the "matching scheme". It has contributions from as 
well as from Cratch = 1 + ci^^ + . . ., namely 



rmatch(^) = -M' - [tT' + ^bQC,]f + 



(14) 



Note that replacing the MS coupling by a non-perturbative one, /match may also be 
defined beyond perturbation theory through eqs. (113181) .^ Another advantage of eq. (fOl) 



^ Clearly the r.h.s. of eq. ( I13> is a function of g^{m\,), i.e. a function of mb/A^. We prefer to write it as 
a function of the ratio of renormalization group invariants, M\, /A^^. 



is that Cps is independent of the arbitrary choice of renormalization scheme for the 
effective operators in the effective theory. Apart from the choice of the QCD coupling, 
the "convergence" of the series eq. ([T?t is dictated by the physics, nothing else. Note 
further that (at leading order in \/M) the conversion function Cps contains the full 
(logarithmic) mass -dependence. The non-perturbative effective theory matrix elements 
are mass independent numbers. Conversion functions such as Cps are universal for all 
(low energy) matrix elements of their associated operator. Thus 

zi(AjWAo(0))QCD [Cps(4)]'(ArW'^r(0))RGi + O(^) (15) 

is a straight forward generalization of eq. (ISJ). 

Analogous expressions for the conversion functions are valid for the time component 
of the axial current replaced by other composite fields, for example the space compo- 
nents of the vector current. Based on the work of lfl2L [l3l [3] and recent efforts their 
perturbative expansion is known including the 3-loop anomalous dimension /match ob- 
tained from the 3-loop anomalous dimension \ 15] and the 2-loop matching function 
lH^flTlflsll . Figured taken from [19.1 . illustrates that the remaining 0(^^(mb)) 



*^match 

errors in Cps seem to be relatively small. 

Although it is generally accepted that 
HQET is an effective theory of QCD in the 
sense that was just described, tests of this 
equivalence are rare and mostly based on phe- 
nomenological analysis of experimental re- 
sults. A pure theory test can be performed 
if QCD including a heavy enough quark can 
be simulated on the lattice at lattice spacings 
which are small enough to be able to take 
the continuum limit. This has recently been 
achieved 1.19.1 and will be summarized below. 
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1 : Cps estimated in perturbation theory. 



2.1. Tests of HQET in a finite volume 

We start with the QCD side of such a test. Lattice 
spacings such that amt, ^ 1 can be reached if one 
puts the theory in a finite volume, xT with L, T 
not too large. We shall use T = L. For various practi- Cb 
cal reasons, so-called Schrodinger functional bound- 
ary conditions are chosen, i.e. Dirichlet in time (at 
xo = 0, T) and periodic in space (id,^. Equivalent 2: The correlation function /a. 
boundary conditions are easily imposed in the effec- 
tive theory [22]. We then form correlation functions of boundary quark fields C, (located 
at xq = 0) and composite fields such as the time component of the axial current in the 
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FIGURE 3. Testing eq. M9\ through numerical simulations in the quenched approximation and for 
L K, 0.2fm L19.I . The physical mass of the b-quark corresponds to z « 5. 



bulk (0 < xo < r), as illustrated in Fig. Eland given for example by 

/a(^o) = -y£((Ai)oWCb(y)75Ci(z)>. (16) 
y,z 

(The current A\ represents the 0(a) -improved version of the axial current for which 
lattice artifacts linear in the lattice spacing are absent.) Another correlation function, 
/i, describes the propagation of a quark-antiquark pair from the xq = boundary to the 
xq = T boundary. For details we refer to fl9] . 

We then take a ratio for which the renormalization factors of the boundary fields 
cancel, 



(a(L)|Ao|5(L)) ^^^^ 



T=L 



|a(L))|||||5(L)) 



|5(L))=e-^^/2|,pB(L)), |a(L))=e-^^/2|(po(L)). (18) 

As shown in the above equations, Ips can be represented as a matrix element of the axial 
current between a normalized state \B{L)) with the quantum numbers of a B-meson and 
\0.{L)) which has vacuum quantum numbers. The time evolution q^^^I^ ensures that 
both of these states are dominated by energy eigenstates with energies around 2/L and 
less. In other words, HQET is applicable if 1/L ^ M (and of course A <^ M). 
One then expects (for fixed LA) 

7ps(L,M)/Cps(M/A)=Xrgi + 0(1/z), z = ML, (19) 

where the l/M corrections are written in the dimensionless variable 1/z and Xrgi is 
defined as Ips but at lowest order in the effective theory and renormalized as in eq. 
Of course such relations are expected after the continuum limit of both sides have been 
taken separately. For the case of Fps(L,M), this is done by the following steps: 

• Fix a value uq for the renormalized coupling g'^{L) (in the Schrodinger functional 
scheme) at vanishing quark mass. In [il9i] uq is chosen such that L ^ 0.2fm. 



• For a given resolution L/a, determine the bare coupling from the condition g^{L) — 
uq. This can easily be done since the relation between bare and renormalized 
coupling is known [23]. 

• Fix the bare quark mass fhq of the heavy quark such that LM = z using the known 
renormalization factor Zm in M = Zmtnq ll23h . 

• Evaluate Yps and repeat for better resolution a/L. 

• Extrapolate to the continuum as shown in Fig. |3l left. 

As can be seen in the figure, the continuum extrapolation becomes more difficult as 
the mass of the heavy quark is increased and 0((aA/)^) discretization errors become 
more and more important. In contrast the continuum extrapolation in the static effective 
theory (Fig.|4l) is much easier (once the renormalization factor relating bare current and 
RGI current is known [6]). After the continuum limit has been taken, the finite mass 
QCD observable yps(L,M) turns smoothly into the prediction from the effective theory 
as illustrated in the r.h.s. figure. Indeed, several such successful tests were performed 
in lHoll . one of them free of the perturbative uncertainty in the conversion function (due 
to reparametrization invariance (2^, ^2^, '26']) and two others with the static (M 0°) 
limit known from the spin symmetry of HQET. For lack of space we do not show more 
examples. We only note that the coefficient of the l/z" terms in naive fits to the finite 
mass results are roughly of order unity. x„„, 
Of course, finite mass lattice QCD results _j ^ 
have been compared to the static limit over the 
years, see for example [27, 3, 4, 28, 29, 30, 31, 
11 S H H S 0] and references therein. 



The new quality of the tests just discussed is 
that the composite fields were renormalized non- 
perturbatively throughout and that, by consider- 
ing a small volume, the continuum limit could be 
taken at large quark masses. 



0.002 0.004 0.006 (a/L)" 

4: Continuum extrapolation of Xrgi fl^ . 



2.2. Beyond the leading order: the need for non-perturbative 

conversion functions 

Both from looking at Fig. Eland from just a naive estimate of A/Mb, one expects that 
the effective theory has to be implemented beyond the leading order in 1/M to reach an 
acceptable precision in this expansion. However, if one wants to do this consistently, i.e. 
one wants to obtain the true coefficients in the expansion, the leading order conversion 
functions such as Cps have to be known non-perturbatively. This general problem in the 
determination of power corrections in QCD is seen by considering the error made in 
eq. S\3[ (or eq. (fTTt ) when the anomalous dimension has been computed at I loops and 
Cmatch at Z — 1 loop Order. The conversion function Cps is then known up to an error 



As mb is made large, this error becomes dominant. Taking a perturbative conversion 
function and adding power corrections to the leading order effective theory is thus to be 
regarded as a phenomenological approach, where one assumes that the coefficient of the 
[g^(mb)]' term is small, such that the (A/Mb)" corrections dominate over a certain mass 
interval. Indeed, returning to our example, Fig.|3lindicates that the power corrections are 
larger than the perturbative ones at 1/z = 0. 1 . . . 0.2. Nevertheless, it remains a fact that a 
theoretically consistent evaluation of power corrections requires a fully non-perturbative 
formulation of the theory including a non-perturbative matching to QCD. 

3. NON-PERTURBATIVE FORMULATION OF HQET 

The discussion in this section summarizes the main points of [8]. We regularize the 
theory on a space-time lattice. In the 1/mb part of the Lagrangian, 

= '£co^'^^^'\x), (21) 

^/') = V7b(-CT-B)v/h, ^i'^ =Whi-i^^)¥h (22) 

the chromo-magnetic field B and the 3-d Laplacian are then discretized in the 

standard way. Details will be irrelevant for our discussion. The coefficients co^^^ are 
functions of the bare coupling go as well as the mass of the heavy quark. They have to 
be determined such as to match the effective theory to QCD. Matching at the classical 
level fixes 

(o[''^ = l/m^ + 0{gl) = 4'^. (23) 

Furthermore, we note that also in eq. ^ a dimension 4 composite field with coefficient 
oc 1/mb -f 0(.?o) added on the r.h.s. when 1/mb corrections are considered. As 

an additional essential ingredient in the formulation of the effective theory we always 

expand the formal weight in the path integral, exp (^Lx — (■^stat(^) + {x) + .. .)^,m 

a power series in 1/mb. The correlation functions are then defined by 

(^) = i|D[(p]^[(p]exp(^-a4£^,tat(x)) x 1 1 -a4£^(i)(x) + . . .| , (24) 

where (p denotes collectively the fields of the theory and the denominator Z insures 
( 1 ) = 1 . The higher order terms in the Lagrangian then appear only as insertions into the 
correlation functions of the static effective theory. The latter is renormalizable by power 
counting and as a result also the effective theory truncated at any finite order in 1/mb 
is renormalizable. With higher dimensional operators in the exponential, as in NRQCD, 
this would not be the case. For the lattice theory renormalizability is important because 
it means that the continuum limit exists and is independent of the details of the lattice 
formulation (universality). 



3.1. Power divergencies 



The coefficients 5m, (O- in eq. @ and eq. (ETT) have a regular expansion in the bare 
coupling gQ. Still, perturbative precision is in general insufficient, since operators of 
higher dimensions mix with those of lower dimension, e.g. 

<='=E^^^r+Ec,^r, c,= '' (25) 
k k ^ 

Since the lattice spacing decreases as a ~ exp( — 1/ (IbogQ)) for small bare gauge cou- 
pling go, a truncation of the perturbative series leaves terms undetermined which diverge 
as the lattice spacing goes to zero. The origin of this problem is the same as the need for 
non-perturbative conversion functions, but the consequence is more drastic due to the 
presence of the hard cutoff in the lattice theory. Without non-perturbative precision for 

5m, COi^\ the continuum limit does not exist. 



3.2. Matching strategy 

The definition of the effective theory is essentially given by eq. (EU, supplemented by 
a definition of the effective composite fields. The only missing piece is a practical strat- 
egy for ascertaining how the parameters in the Lagrangian and in the effective fields can 
be determined beyond perturbation theory. At a given order n in the 1/mb-expansion, we 
denote the parameters in the effective theory hy Cj^,k= I, ... ,Nn. Observables, e.g. renor- 
malized correlation functions or energies are denoted by ^»^^^^(L,M) (4>^*^'^(L,M)) in 
the effective theory (in QCD), with the argument M referring to the mass of the heavy 
quark and L the linear extent of the finite volume. It is then sufficient to impose 

4>«Q^T(Lo,M)=4>QC°(Lo,M), k=l,...,N„. (26) 

to determine all parameters {ck , A: = 1 , . . . , A^„} in the effective theory. Observables used 
originally to fix {c^.k = l,...,A/f}, the parameters of QCD, may be amongst these 
^QCD rpj^g matching conditions, eq. (l26t . define the set {q} for any value of the lattice 
spacing (or equivalently bare coupling). Here, a typical choice is Lq ^ 0.2...0.4fm, 
since in such a volume the r.h.s. of the equation can be evaluated, see Sect. 12.11 In 
practice, the parameters of the effective theory are then determined at rather small lattice 
spacings in a range of a = 0.01 fm to a = 0.04fm. Large volumes as they are needed to 
compute the physical mass spectrum or matrix elements then require very large lattices 
(L/a > 50). A further step is needed to bridge the gap to practicable lattice spacings. A 
well-defined procedure is as follows: First we assume that all observables 4>"'^^^(L,M) 
have been made dimensionle ss by multiplication with appropriate powers of L. Next, we 
define step scaling functions ll38l] . F^, by 

4>«Q^^(5L,M) = F,({4>«Q^T(L,M) ,7 = 1,... , A^„}) , = 1, . . . , , (27) 



where typically one uses scale changes of = 2. These dimensionless functions describe 
the change of the complete set of observables {^^*^^^} under a scaling of L ^ sL. 
In order to compute them one selects a lattice with a certain resolution a/L. The 

specification of ^>^^^^(L,M), j = l,...,Nn, then fixes all (bare) parameters of the 
theory. The l.h.s. of eq. dTTl is now computed, keeping the bare parameters fixed while 
changing L/a — L' /a = sL/a. The values for the continuum can then be be reached 
by extrapolating the resulting lattice numbers to a/L ^ 0. 

After repeating this step two or three times with s = 2, lattice spacings appropriate for 
infinite volume computations will be reached. 



3.3. Example: the mass of the b-quark at lowest order 

For illustration purposes we consider a simple example here, the computation of 
the b-quark mass, starting from the observed B-meson mass. Already at the lowest 
order in 1/mb the mixing of operators of different dimensions is relevant in this case: 
"y7jjDoV4i mixes with V/^V^i- Hence dm = {cigQ + C2gQ + ...)/ a, the coefficient of xfr^xir^ 
in the Lagrangian eq. ©, has to be determined (or eliminated) non-perturbatively. In 
eq. (l26b we have n = 0, Nq = Nf+ I and we omit the discussion of the choices for 4>fc, 
k = 1, . . . , A^f which fix the bare light quark masses and coupling (both in QCD and 
HQET). Qbviously any energy in the b-sector will do to fix 5m. We choosefl8l l39ll 

r(L,M) = ^ln[/A(.^o-a)//A(xo + a)] (V^ fixed) (28) 

and we require eq. (E^ where for ^ = A^f + 1 we identify 

^f+,mM) = LT{LM) , 4>«Qf (L,M) = L(rstat(L) +m) , (29) 

Here m represents the quark mass that was removed from all energies in defining the 
effective theory such that the m^ oo limit exists and Fstat refers to eq. (ESI at the lowest 
order in l/M. The relevant part of eq. (l27t can then be written in the simple form, 

4>«Q7(2L,M) = 2^^f^^{L,M)+am{f{L)), (30) 

OmifiL)) = 2L[rstat(2L)-rstat(L)]. (31) 

In Om the divergent 5m (as well as the mass shift m) cancel. It is independent of the 
mass. We now see immediately that 

me = fitat - rstat(L2)^ + rstat(L2) - rstat(Lo)^ + r(Lo,Mb) +Q(AVMb) . (32) 

a^o in HQET a^o in HQET a^o in QCD 

Here, Lgtat = hniL^cx, Fstat (L) is the infinite volume energy of a B-meson in static ap- 
proximation. It is often called the static binding energy. The whole strategy is illustrated 



In practice F is replaced by the spin averaged energy is Hsoll . 



in Fig. 13 As indicated in eq. (B^ . the continuum limit can be taken in each individual 
step; a numerical example is shown in Fig.0 
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FIGURE 5. Connecting experimental observables to renormalized HQET. Fstat is a renormalized 
quantity in HQET and am{g'^{L)) connects Fstat (i) and Fstat (2L). 

After obtaining all pieces in eq. (l32b . the equation is numerically solved for Zb = M^Lq. 
Since also the size of Lq in units of vq ^ 0.5fm lEol] is known, one can quote (remember 
mb is in the MS scheme at scale m\,) 



roMb= 16.12(24)(15) ^ Mb = 6.36(10)(6)GeV, mb = 4.12(7)(4) GeV. (33) 

We emphasize that this result is in the quenched approximation but includes the lowest 
non-trivial order in 1/mb. An estimate of the associated 0(A^/Mb) uncertainty is not 
included in the errors shown. Our discussion mainly serves to illustrate the potential 
of the approach in an example where the power divergent mixing needed to be solved 
non-perturbatively. 



4. PERSPECTIVES 



Non-perturbative HQET at the leading 
order in 1/m has reached a satisfac- 
tory status, with the b-quark mass 
and the Bs-meson decay constant [6, 41] 
known in the continuum limit of the 
quenched approximation. Their preci- 
sion can and will still be improved. Ap- 
plying these methods to the theory with 
dynamical fermions is straight forward; 
"only" the usual problems with the light 
quarks have to be solved. By themselves 
such lowest order results are not ex- 
pected to have an interesting precision 
for phenomenological applications, but 

certainly they can constrain the large mass behavior computed with other methods 
121 J„ A 2a, 22, ia IL i2, i3, i4, i6, . 

More interestingly, also 1/mb corrections can, in principle, be computed in the effec- 
tive theory. Here, details of the necessary numerical steps have not yet been implemented 
but the very first tests have been encouraging [42] . Another relevant technical advance 
has been the realization that a change of the regularization details allows to achieve much 
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better statistical errors in HQET, while keeping the discretization errors small (4?]. In 
summary, we believe that all ingredients exist which are needed to apply HQET beyond 
the leading order in 1 /mi,. 
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